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ABSTRACT 

^ ! A double-peak spectral-energy-density of 7-rays, similar to that observed in 

Q ! blazars, is expected in 7-ray bursts (GRBs) produced in supernova (SN) explo- 

^ i sions. The second peak, which is formed by inverse Compton scattering of the 

O ! ambient SN light by cosmic ray electrons accelerated by the jets from the SN 

O I explosion, has a much higher peak-energy than the first ordinary peak. However, 

^ I in X-ray flashes (XRFs), which in the cannonball (CB) model are normal GRBs 

I viewed farther off axis, the first peak-energy is shifted to the soft X-ray region 

^ I while the second peak-energy moves to the MeV range. In far-off-axis GRBs, 

\0 [ such as GRBs 980425 and 031203, the second peak may have been confused 

^ I with the normal GRB peak. In most GRBs, which have been observed so far, 

^ I the 7 ray detectors ran out of statistics far below the second peak. However, 

Tjj- I in bright GRBs, the two peaks may be resolved by simultaneous measurements 

O ■ with SWIFT and GLAST. 

^ : 
o : 

^ , 1. Introduction 

■ 

• There is mounting evidence that long duration 7-ray bursts (GRBs) and X-ray flashes 

■ (XRFs) are both produced by highly relativistic and narrowly collimated jets ejected in 

' supernova (SN) explosions akin to SN1998bw (see, e.g. Dar 2004 and references therein). 
■ - - ' Low-luminosity GRBs and XRFs seem to be ordinary GRBs whose 7-rays are much softer 
and their luminosity is much lower because they are viewed from angles relative to the jet 
which are a few times larger than the viewing angles of ordinary GRBs (see, e.g. Dar & 
De Riijula 2000,2003; Dado, Dar & De Rujula 2002,2004 and references therein). However, 
recently, the European Space Agency satelhte INTEGRAL discovered GRB 031203 (Gotz 
et al. 2004) at a redshift z = 0.1055 (Prochaska et al. 2004), the second nearest GRB after 
GRB 980425 at z = 0.0085 (Galama et al. 1998), which hke GRB 980425 had extremely 
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low lunminosity, but, contrary to the expectation, its 7-rays were not much softer than 
those of ordinary GRBs (Sazonov et aL 2004). GRB 031203 was produced in an SN very 
similar to SN1998bw which produced GRB 980425 (e.g. Malesani et al. 2004 and references 
therein). Consequently, the INTEGRAL observations were interpreted as evidence of a 
different class of GRBs, which are intrinsically very faint and include GRB 980425 and GRB 
031203 (Sazonov et al. 2004; Soderberg et al. 2004; Woosley et al. 2004) and are produced 
in SNe akin to SN1998bw. 

The spectrum of the soft 7-rays emitted in GRB 031203, which was measured with 
INTEGRAL, is well described by dn^/dE ~ £;-i-63±o.o6 between 20 and 400 keV (Sazonov et 
al. 2004). But, its soft X-ray fiucncc in the 0.2-10 keV range, - 2.6 ±1.3 x lO''^ erg s'^ , 
which was inferred (Watson et al. 2004) from its X-ray dust-scattcrcd echo measured by 
XMM-Newton (Vaughan et al. 2004), is far above the INTEGRAL spectrum extrapolated 
to the soft X-ray region. It has been suggested that the soft X-ray fluence may have been 
overestimated (Prochaska et al. 2004; Sazonov et al. 2004). 

However, in this letter we suggest an alternative interpretation of the INTEGRAL and 
XMM-Newton observations, which is based on the cannonball (CB) model of GRBs and 
XRFs (e.g. Dar & De Riijula 2003; Dado et al. 2004) and which removes the conflicts between 
the off-axis interpretation of the unusually low luminosity of GRB 031203 and its high peak- 
energy, and between its X-ray fluence extrapolated from the INTEGRAL measurements 
(Sazonov et al. 2004) and that inferred from the XMM-Newton measurements (Vaughan et 
al. 2004). We also demonstrate that this interpretation can explain other puzzhng GRB 
observations. 

The CB model assumes that in core-collapse SN, after the collapse of the stellar core 
into a neutron star or a black hole, an accretion disk or a torus is produced around the newly 
formed compact object, either by stellar material originally close to the surface of the im- 
ploding core and left behind by the explosion-generating outgoing shock, or by more distant 
stellar matter falling back after its passage (De Riijula 1987). A highly relativistic plasmoid 
(CB) of ordinary matter is emitted along the rotation axis when part of the accretion disk 
falls abruptly onto the compact object, as observed in microquasars (e.g. Mirabel & Rodrigez 
1999; Rodriguez & Mirabel 1999 and references therein). A CB contains a thermal plasma 
with a power- law tail of knocked-on electrons from collisions with particles of the interstellar 
medium (ISM) and swept-in ISM electrons, which are Fermi-accelerated and cool quickly 
by synchrotron emission from the strong equipartition magnetic field in the CB to a steady 
state distribution, dn^/dE ~ with p ^ 2.2, in the CB rest frame. As the jet of 

CBs coasts through the "ambient light" permeating the surroimdings of the parent SN, the 
electrons enclosed in the CBs Compton up-scatter photons to energies which, close to the 
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CBs direction of motion, correspond to the 7-rays of an ordinary GRB and less close to it, 
to the X-rays of an XRF (Dar & De Riijula 2003). 

A CB also produces a narrow conical jet of high energy cosmic-ray electrons (and nuclei) 
along its motion (e.g., Dar & De Riijula 2003) . These electrons are of two origins: swept-in 
ISM electrons which were Fermi-accelerated in the CB and escaped out into the ISM and ISM 
electrons which were scattered elastically by the CB. These cosmic-ray electrons produce a 
second peak in the spectral energy flux of GRBs and XRFs at a much higher energy, by 
inverse Compton scattering (ICS) of "ambient light" permeating the surroundings of the 
parent SN. 

Thus, the internal CB electrons and the external cosmic ray electrons produce, by ICS 
of ambient hght, a spectral-energy-density {E'^ dn/dE) of 7-rays with two peaks, like that 
observed in blazars (e.g., Padovani & Giommi 1995; Wehrle et al. 1998). Normal GRBs 
have their first peak-energy usually around a fraction of an MeV and then their second peak- 
energy at a much higher energy. Gamma-ray detectors on board satellites usually ran out of 
sensitivity/statistics well below the second peak-energy. However, high-energy photons with 
a flux much larger than that expected from the extrapolated decline of the first peak, have 
been discovered in a few cases of very bright GRBs with instruments on board CGRO (e.g. 
Hurley et al. 1994; Dingus et al. 1995; Gonzalez et al. 2003). In this letter we demonstrate 
that such a high energy component is well described by the CB model. Moreover, in very 
low-luminosity GRBs/XRFs, which are far off-axis GRBs, the first peak is in the keV range 
while the second peak is in the MeV range, i.e., the spacing between the peaks is relatively 
small and both peaks, or a significant fraction of both, may fall within the detection range 
of the ordinary XRF and GRB detectors. A score of low-luminosity XRFs were detected 
by CGRO, BeppoSAX and HETE 11. But, in the case of HETE II, its FREGATE detector 
could go only to 400 keV. BATSE on board CGRO and GRBM on board BeppoSAX usually 
ran out of statistics/sensitivity well below the second peak, although a significant emission 
in excess of that expected from a simple power-law extrapolation of the low energy spectrum 
to higher energies was detected in some XRFs (see e.g., Frontera et al. 2004a,b). In this 
letter we first demonstrate that the peak-energy and the isotropic radiation energy emitted 
in the first peak satisfy well the simple correlation (1 -|- ^) £'p ex , predicted by the 

CB model. GRBs 980425 and 031203, which appear to be exceptions, could have been 
far-off-axis GRBs with a relatively small spacing between their keV and MeV peaks, and 
therefore appeared as GRBs with a '"normal" peak-energy. We predict that simultaneous 
measurements with SWIFT and GLAST of the spectral energy fiux of ordinary GRBs will 
resolve it into two peaks which shift to smaller energies during the GRB pulse. 
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2. The first peak 

Let 7 be the Lorentz factor of a CB and S — 1/7(1-/3 cos^) be its Doppler factor 
when viewed from an angle 9 relative to its motion. In the CB model (Dar & De Riijula 
2003), the observed peak-energy of 7-rays produced at a redshift z by ICS of thin thermal 
bremsstrahlung light around an SN, dn^/dE ~ exp{—E/T)/E with a typical energy e-y ~ 
T ~ 1 eV^ , is given by, 

Under the assumption of isotropic emission in the CB rest frame, Doppler boosting and 
relativistic beaming yield a 7-ray fluence of a GRB pulse, which is proportional to 5^ , 

where E'^ is the total 7-ray energy emitted in the CB rest frame. Consequently, El^°, the 
inferred 'GRB isotropic 7-ray energy' in a GRB pulse, under the assumption of isotropic 
emission in the observer frame, is 



£;^°=i!^F^f«5^£;;. (3) 



In the CB model, the predicted GRB spectrum from ICS of ambient light with a thin 
thermal bremsstrahlung spectrum, by the electrons inside the CB, is. 




where Te// — ^5T/{1 -\- z) and 6 is a dimensionless constant. The values of a and (5 
may deviate from their indicated values, because the ambient radiation may deviate from 
a thin thermal bremsstrahlung, and the power-law index of the accelerated and knocked-on 
electrons after cooling may be larger than p + 1 = 3.2 and increase with time. Also the 
temperature of the ambient light seen by the CB is decreasing with distance, approximately 
as 

T,ff{t) ~ Te;^(0) {1 - exp[-{t,/tf]}, (5) 

where to is a constant. As was shown in Dar & De Rujula 2003 and in Dado et al. 2004, 
Eq. (4) is practically indistinguishable from the phenomenological Band function (Band et 
al. 1993) and it is in good agreement with the measured spectrum of the photons in the first 
peak of the spectral-energy-density of ordinary GRBs and XRFs. 
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If core collapse SNe and their environments were all identical, and if their ejected CBs 
were also universal in number, mass, Lorentz factor and velocity of expansion, all differences 
between GRBs would depend only on the observer's position, determined by z and the angle 
of observation, Q. For a distribution of Lorentz factors that is narrowly peaked around 
7 ~ 10^, the ^-dependence is in practice the dependence on 5, the Doppler factor. Hence 
Eqs. (1),(3) yield the correlation (Dar and De Rujula 2000; 2003), 



In Fig. 1 we plot (1 + z)Ep as function of E\^° for GRBs/XRFs of known rcdshift, as 
compiled by Amati (2004) and by Ghiralanda et al. (2004), and the best fitted power-law, 
{\^ z)Ep = a [Ei^""]'^ . The best fitted power-law index obtained with the CERN standard 
program MINUIT, a = 0.335 ± 0.06 and a = 0.345 ± 0.07, for the Amati (2004) and 
Ghiralanda et al. (2004) comilations, respectively, is in agreement with Eq. (6) but not with 
the original relation, (l + z) Ep '-^ a [£;«'Soj0.52±o.06 proposed by Amati et al. (2002), which, as 
far as we know, has no theoretical derivation. We also plot in Fig. 1 parallel lines between 
which the observed values are expected to be spread because of the spread in the 'standard 
candle' properties. The correlation predicted by the CB model is well satisfied except for 
GRB 980425 and GRB 031203. 



In the CB model, the ISM in front of the CBs is ionized by the beamed radiation from 
the highly relativistic CBs. The turbulent magnetic fields in the CBs accelerate the ionized 
ISM particles, which they gather on their path, to an initial distribution, drie/dE ~ E~p 
with p ~ 2.2 in the CB rest frame. In a steady state situation, the electrons, which are 
trapped in the CB by its internal magnetic fields and cool quickly by synchrotron emission, 
reach a distribution, dUg/dE ~ E~'^'^^\ while the electrons which escape the CB must have 
the hard "injection spectrum" , due/ dE ~ E'P . Their cooling time in the ISM is much longer 
than that of the electrons which are trapped magnetically inside the CBs, because the ISM 
magnetic field is smaller by many orders of magnitude than that inside the CBs. 

In the CB rest frame, the Lorentz boosted ambient light undergoes Compton scattering 
from these two distributions and produces the first peak with the "Band function" shape 
and a hard tail at higher energies with practically a "time- independent" spectral index 1.6, 



(i + z)Epo^[Ei;"]y\ 



(6) 



3. The second peeik; 



dn. 



-1.6 



(7) 



dE 



The ISM electrons which are scattered elastically by the highly relativistic CBs in the 
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direction of the observer have an approximate lab energy, Eg ^ 'y 6 rrie ■ Because of their very 
large Lorentz factor, Fg ~ 7 5 , the ambient photons which they scatter have much higher 
energies than those of the photons scattered by "cold" electrons in the CBs, and they are 
narrowly beamed along the electrons 'direction of motion. Hence, in the Thomson regime 
they have approximately the thin thermal bremsstrahlung distribution of the ambient light 
boosted by ~ 4 7^ 5^/3 : 

^^o.Ee-'^/'^^ff, (8) 

where T^ff = (4/3) 7^ 5^ T/ (l + z) . For E -C T^fj , the low energy side of the spectral energy 
density of this "second peak" has the simple shape, E'^ dN^/dE oc E . 

During the GRB phase, the radius of a CB increases linearly with the distance x from 
the explosion site while the density of the ambient light decreases like l/x^ . Because of the 
slow cooling rate in the ISM, the accumulated number of high energy electrons in the narrow 
beam produced by the highly relativistic CB is proportional to the ISM mass swept up by 
the CB. Thus, for a constant density-profile, the magnitude of the spectral energy density, 
which is proportional to the product of the swept-up mass and the density of the ambient 
light, increases linearly with distance, whereas it decreases like l/x for a density-profile of a 
stellar wind which blows constantly and produces a density ng oc . During the GRB the 
deceleration of the CBs is negligible and the observer time is proportional to the distance. 
Thus, for a constant density-profile the spectral energy density increases with time like t , 
whereas it decreases like 1/t for a "windy" density-profile. 

In the Thomson regime, the peak-energies are given approximately by Ep Tgff/{l+z) . 
Consequently, Ep[l] and Ep[2] , the peak energies of the first and second peak, respectively, 
are related through 

^j2l.Ii±M. (9) 

For ordinary GRBs with Ep[l] ~ 200 keV and z ~ 1 , Eq. (9) yields Ep[2] ~ lOOGeV , 
whereas for a very dim XRF with Ep[l] ~ IkeV and z « 1 , it yields, Ep[2] ~ IMeV . 
During each CRB pulse, both peak-energies, being proportional to T{t) , decrease during a 
GRB pulse like, Teff{0) [1 - exp[-{to/tf]] . 



4. Comparison with observations 

There are very few cases where GRBs were measured over a very broad energy range 
and allow a comparison between the CB model predictions and the observations. Here we 
shall limit the comparison to two GRBs of particular interest. 
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4.1. 



GRB 941017 



Gonzalez et al. (2003) reported the discovery a high energy spectral component in 
GRB 940117 with an energy flux density, E^dn/dE - E, peak energy Ep[2] > 200 MeV 
and a fluence > 3 times that of the the first normal GRB component. Their results are 
shown in Fig. 2 borrowed from Gonzalez et al. 2003. The theoretical lines are best fitted 
Band functions plus a power-law contribution (Gonzalez et al. 2003). But, these lines are 
indistinguishable from the bes-fitted CB model double-peak spectral-energy-density which 
is a sum of E'^ dn/dE with dn/dE as given by Eq. (4) and of a second component which is 
described by Eq. (8) with E < Te//[2] . 



GRB 031203 (Gotz et al. 2003) was produced in an SN explosion similar to SN1998bw 
(Bersier et al. 2004; Cobb et al. 2004; Malesani et al. 2004; Thomsen et al. 2004; Gal- Yam 
et al. 2004) and had an unusually low inferred isotropic luminosity. Its soft X-ray emission, 
which was inferred from modehng its measured dust scattered echo by XMM- Newton, yielded 
(Watson et al. 2004; Vaughan et al. 2004)), ~ 2.6± 1.3 x 10"^ erg cm-^ , in the 0.7-5 keV 
range, far above the extrapolation of the INTEGRAL spectrum. This can be seen from the 
upper panel of Fig. 3 where we demonstrate a CB model fit to the spectral energy density 
of GRB 031203 which consists of two terms, E"^ dn/dE with dn/dE as given by Eq. (4) for a 
low-energy peak with Ep ^ T^ff ^ 3 keV and b = 10^^ , and a hard-tail contribution as given 
by Eq. (7). The normalization of both components were adjusted to fit the observational 
data. In the lower panel of Fig. 3b we compare the photon spectral index of the hard tail 
component, (p+ l)/2 Ri 1.6, as predicted by the CB model, Eq. (7), and its measured value 
by INTEGRAL (Sazonov et al. 2004). 



In the CB model, the internal population of electrons produce the ordinary peak in 
the spectral-energy-density of GRBs/XRFs and the external population of high energy (cos- 
mic ray) electrons, which are accelerated by the CBs, produce a second GRB peak at a 
much higher energy through inverse Compton scattering of ambient light permeating the 
surroundings of the parent SN. We have shown that for GRBs with known redshift, usually 
the peak-energy and the isotropic radiation energy of the first, ordinary peak satisfy well 
the simple correlation {1 + z) Ep oc [E!^'"]^/^ , predicted by the CB model. However, GRBs 



4.2. 



GRB 031203 



5. 



Conclusions 
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980425 and 031203 have peak-energies much larger than those expected in the CB model 
from their very low luminosities. But, in very low luminosity GRBs and XRFs which, in 
the CB model, are ordinary GRBs viewed far off axis, the second peak moves down to the 
MeV region and can dominate the soft 7-ray region. We have shown that GRBs 031203 and 
GRB 980425 could have been such far-off-axis GRBs with their first peak in the soft X-ray 
region but with a contribution from their second broad peak which dominates their X-ray 
and 7-ray emission, as shown here exphcitely for GRB 031203. We have also shown that 
the high energy component which was discovered in GRB 941017 is well explained by the 
second peak. Finally, wc predict that simultaneous measurements with SWIFT and GLAST 
will resolve the spectral energy flux of ordinary GRBs/XRFs into a double peak spectrum 
whose peaks move to lower energies during the GRB pulse. 

Acknowledgement: Useful comments by Avishai Gal- Yam, Kevin Hurley and Simon 
Vaughan are gratefully acknowledged. This researcch was supported in part by the Asher 
Fund for Space Research at the Technion. 

REFERENCES 

Amati, L., et al. 2002, A&A, 390, 81 

Amati, L. 2004, astro-ph/0405318 

Band, D., et al. 1993, ApJ, 413, 281 

Bersier, D., et al. 2004, GGN Girc. 2544 

Gobb, B. E., et al. 2004, astro-ph/0403510 

Dado, S., Dar A. & De Riijula, A. 2002, A&A, 388, 1079 

Dado S., Dar A. & De Riijula A. 2004, A&A, 422, 381 

Dar, A. 2004, astro-ph/0405386 

Dar, A. & De Riijula, A. 2000, astro-ph/0012227 

Dar, A. & De Riijula, A. 2003, astro-ph/0308248 (Phys. Rep. in press) 

De Riijula, A. 1987, Phys. Lett. 193, 514 

Dingus, B. L., et al. 1995, Ap&SS, 231, 187 

Frontera, F., et al. 2004b, astro-ph/astro-ph/0406575 



-9- 

Frontera, F., et al. 2004, astro-ph/0408436 
Galama, T. J., et al. 1998, Nature, 395, 670 
Gal- Yam, A., et al. 2004, ApJ, 609, L59 
Gonzalez, M. M., et al. 2003, Nature, 424, 749, 
Gotz, D., et al. 2003, GCN 2459 

Ghirlanda, G., Ghisellini, G. & Lazzati, D, 2004, astro-pli/0405602 
Hurley, K., et al. 1994, Nature, 372, 652 
Malesani, D., et al. 2004, astro-ph/0405449 
Padovani, P. & Giommi P. 1995, ApJ, 444, 567 
Prochaska, J. X., et al. 2004, ApJ, 611, 200 

Sazonov, S. U., Lutovinov, A. A., & Sunyaev, R. A. Nature, 2004, 430, 646 

Soderberg, A. M., et al. 2004, Nature, 430, 648 

Thomsen, B., et al. 2004, astro-ph/0403451 

Vaughan, S., et al. 2004, ApJ, 603, L5 

Watson, D., et al. 2004, ApJ, 605, LlOl 

Wehrle, et al. 1998, ApJ, 497, 178 

Woosley, S., 2004, Nature, 430, 623 

Yamazaki, R., Yonetoku, D. & Nakamura, T. 2003, ApJ, 594, L79 



This preprint was prepared with the AAS IMJ^pC macros v5.2. 



10"^ 10"^ 1 10 10^ lo' 

Eiso (10^' erg) 




Fig. 1. — The observed peak-energy corrected for redshift, {1 + z) Ep as function of the 
inferred isotropic radiation energy emitted at redshift z for GRBs/XRFs of known redshift 
as compiled by (a) Amati (2004) and by (b) Ghiralanda et al. (2004). The thick hnes are the 
best fitted power-laws: A £;0-334±o.06 ^^^^ ^ ^o.345±o.07 ^ respectively. The dotted lines border 
the estimated spread (a factor of ~ 4) in the isotropic radiation energy due to the spread in 
the number of CBs and their physical properties in GRBs. The point due to GRB 980425 
(as inferred by Yamazaki et al. 2003) is far above the theoretical expectation. 




Fig. 2. — The spectral energy density of GRB 941017 obtained by Gonzalez et al. (2003) 
with BATSE-LAD (crosses) and EGRET-TASC (filled circles) on board CGRO. The CB 
model fit for the first peak is indistinguishable from the "Band function" fit (dashed line) 
by Gonzalez et al. (2003). The rising side of the second peak is well described by the CB 
model prediction as given by Eq. (8). The fits are shown for the time intervals (-18)- 14s (a); 
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Fig. 3. — Top: Comparison between the observed spectral energy density of GRB 031203 
and a CB model fit given by Eq. (4) with Ep f» T^ff ^ 3 keV and b = 10^^ , and a "hard tail" 
as given by Eq. (7). The X-ray fiuence is that inferred by Vaughan et al. (2004) from the 
dust echo observed with XMM-Newton. The INTEGRAL measurements are those reported 
by Sazonov et al. (2004). Bottom: Comparison between the the photon spectral index 
as obtained from single power-law fit to the photon spectrum measured with INTEGRAL 
(Sazonov et al. 2004) and the CB model prediction for the hard-tail component in GRBs 
generated by cosmic ray electrons which were accelerated in the CBs. 



